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Abstract 

Neutrinos coupled to an underlying scalar field in the scenario for unification of mass varying 
dark matter and cosmon-Zz&e dark energy is examined. In the presence of a tiny component of mass 
varying neutrinos, the conditions for the present cosmic acceleration and for the stability issue are 
reproduced. It is assumed that sterile neutrinos behave like mass varying dark matter coupled 
to mass varying active neutrinos through the seesaw mechanism, in a kind of mixed dark matter 
sector. The crucial point is that the dark matter mass may also exhibit a dynamical behavior 
driven by the scalar field. The scalar field mediates the nontrivial coupling between the mixed 
dark matter and the dark energy responsible for the accelerated expansion of the universe. The 
equation of state of perturbations reproduce the generalized Chaplygin gas (GCG) cosmology so 
that all the effective results from the GCG paradigm are maintained, being perturbatively modified 
by neutrinos. 
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I. INTRODUCTION 



The continuous search for the physical mechanism which sets the cosmic acceleration 
of the universe and the corresponding conditions for stability has stimulated interesting 
and sometimes fascinating discussions on cosmological models [1-51. The dynamical mass 

na 

attributed to neutrinos or dark matter [y, |7J can, for instance, regulate the time evolution 
of the dynamical dark energy providing the setup of the cosmic acceleration followed by 
the cosmological stability In this context, the coupling of mass varying dark matter with 
neutrinos yields interesting relations between the present mass of neutrinos and the dark 
energy equation of state. 

In a previous issue js], it was demonstrated that an effective generalized Chaplygin gas 
(GCG) scenario js-ll] can be reproduced in terms of a dynamical dark energy component 
4> with equation of state given by p<j,{$) = —p^) and a cold dark matter (CDM) component 
with a dynamical mass driven by the scalar field 0. Dark matter is, most often, not con- 
sidered in the mass varying neutrino (MaVaN) models. The treatment of dark energy and 
dark matter in the GCG unified scheme naturally offers this possibility. Identifying sterile 
neutrinos as dark matter coupled with dark energy provides the conditions to implement 
such unified picture in the MaVaN formulation since active and sterile neutrino states are 
connected through the seesaw mechanism for mass generation. The constraints imposed by 
the seesaw mechanism allows one to establish an analytical connection to the GCG in terms 
of a real scalar field. The dynamics of the coupled fluid composed by neutrinos, dark matter 
and dark energy is driven by one single degree of freedom, the scalar field, 0(a). 

The simplest realization of MaVaN mechanisms 15Ml9| consists in writing down an ef- 
fective potential which, in addition to a scalar field dependent term, contains a term related 
to the neutrino energy density. It results in the so-called adiabatic condition, which some- 
;imes leads to the stationary regime for the scalar field in respect with an effective potential 
171 . |2CH22| . One indeed expect a tiny contribution from cosmological neutrinos to the en- 
ergy density of the universe. MaVaN scenarios essentially predict massless neutrinos until 
recent times. When their mass eventually grows close to its present value, they form a 
non-relativistic (NR) fluid and the interaction with the scalar field stops its evolution. The 
relic particle mass is generated from the vacuum expectation value of the scalar field and 
becomes linked to its dynamics: m(4>). It is presumed that the neutrino mass has its origin 
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on the vacuum expectation value (VEV) of the scalar field and its behavior is governed by 
the dependence of the scalar field on the scale factor. 

In fact, it is well-known that the active neutrino masses are tiny as compared to the 
masses of the charged fermions. This can be understood through the symmetry of the stan- 
dard model (SM) of electroweak interactions. It involves only left-handed neutrinos such that 
no renormalizable mass term for the neutrinos is compatible with the SM gauge symmetry 
SU (2)l<S>U(1)y- Once one has assumed that baryon number and lepton number is conserved 
for renormalizable interactions, neutrino masses can only arise from an effective dimension 
five operator. It involves two powers of the vacuum expectation value of the Higgs doublet. 
They are suppressed by the inverse power of a large mass scale M of a sterile right-handed 
Majorana neutrino, since it has hypercharge null in the SM. This super-massive Majorana 
neutrino should be characteristic for lepton number violating effects within possible exten- 
sions beyond the SM. At our approach, the mass scale M has its dynamical behavior driven 
by 0. In some other words, the sterile neutrino mass is characteristic of the aforementioned 
mass varying dark matter which indirectly results in MaVaN's, i. e. active neutrinos with 
mass computed through seesaw mechanism. 

In section II, we report about the main properties of a unified treatment of dark matter 
and dark energy prescribed by the mass varying mechanism. One sees that model dependent 
choices of dynamical masses of dark matter allows for reproducing the conditions for the 
present cosmic acceleration in an effective GCG scenario. The stability conditions resulted 
from a positive squared speed of sound, c 2 s , is recovered. In section III, we discuss the neutrino 
mass generation mechanism in the context of the GCG model. Following the simplest 
formulation of the seesaw mechanism, Dirac neutrino masses with analytical dependencies 
on 0, 1 and _1 are considered. In section IV, we discuss the conditions for stability and the 
perturbative modifications on the accelerated expansion of the universe in the framework 
here proposed. One can state the conditions for reproducing the GCG scenario. We draw 
our conclusions in section V. 

II. MASS VARYING DARK MATTER AND THE GCG 

To understand how the mass varying mechanism takes place for different particle species, 
it is convenient to describe the relevant physical variables as functionals of a statistical 
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distribution f(q). This counts the number of particles in a given region around a point of 
the phase space defined by the conjugate coordinates: momentum, (3, and position, x. The 
statistical distribution f(q) can be defined in terms of the comoving momentum, q, and for 
the case where f(q) is a Fermi-Dirac distribution, it can be written as 

f( q ) = {exp[(q-q F )/T } + l}- 1 , (1) 

where T is the relic particle background temperature at present. 

In the flat FRW scenario, the corresponding particle density, energy density and pressure 
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can thus be depicted from the Einstein's energy-momentum tensor 

i r 

7r 2 a 3 J 

PMiaA) = / d ^ 2 (q 2 + M 2 Wa 2 ) l/2 f( q ), (2) 

7r 2 a 4 ' 
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PM(aA) = ^r- 4 I dqq 4 (q 2 + M 2 (<f>) a 2 y 1,2 f{q). 

07T a 







where a is the scale factor (cosmological radius) for the flat FRW universe, for which the 
metrics is given by ds 2 = dt 2 — a 2 (t)S ij dx i dx j . 

From the dependence of p on a one can obtain the energy-momentum conservation equa- 
tion, 

p M + 3H(p M + Vm) ~ i>-^-^ = °> ( 3 ) 

that translates the dependence of M on into a dynamical behavior, and where H = a/a 
is the expansion rate of the universe and the overdot denotes differentiation with respect to 
time (' = d/dt). Simple mathematical manipulations allow one to easily demonstrate that 

M ( a )^77T = - 3 PmM), (4) 

OlVl (a) 

from which, one can see that the strength of the coupling between relic particles and the 
scalar field is suppressed by the relativistic increase of pressure (p ~ 3p). As long as particles 
become ultra-relativistic (T(a) = T /a oc q/a » M(4>(a))) the matter fluid and the scalar 
field fluid tend to decouple and evolve adiabatically separated. 

Adding Eq. ()3]) to the equation of motion for a cosmon field <ft like 

^Min^M dUW n,^ dM dpM 
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rewritten as 



+ 3H(o 6 + 



■dMdp M 
+ (P^-^r^ = 0, 



(6) 



d<P dM 

results in the equation of energy conservation for a unified fluid with a dark energy compo- 
nent and a mass varying dark matter component, 



p + 3H(p + p) = 0. 



(7) 



for which we identify p = p^ + Pm and p = p^ + Pm- 

From this point, it is suggested that such a unified fluid corresponds to an effective 
description of a GCG universe. The GCG model is characterized by an exotic equation of 



state 
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12j given by 



p = -A s p 



(8) 



which can be obtained from a generalized Born-Infeld action Inserting the above equa- 
tion of state into the unperturbed energy conservation Eq. ([7]) and following a straightforward 
integration [llj], one obtains 



P = Po 



A.) 



i/(i+«) 



and 



P 



-A s p 



A, 



a s(i+a) 
(1-^) 



*/(l+a) 



(9) 



(10) 



a s(l+a) 

One of the most striking features of the GCG fluid is that its energy density interpolates 
between a dust dominated phase, p oc a~ 3 , in the past, and a de-Sitter phase, p = —p, at 
late times. This property makes the GCG model an interesting candidate for the unification 
of dark matter and dark energy. Notice that for A s = 0, GCG behaves always as matter 
whereas for A s — 1, it behaves always as a cosmo logical constant. Hence to use it as a 
unified candidate for dark matter and dark energy one has to exclude these two possibilities 
so that A s must lie in the range < A s < 1. Furthermore, this evolution is controlled by 
the model parameter a. 

Assuming the canonical parametrization of p and p in terms of a scalar field 4>, 

1 •„ 



P 
P 



2 : 2 + 



(11) 
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and following Ref. [13J, one can obtain through Eq. (l9] fT0]) the effective dependence of on 
a implicitly given by 



2 (a) 



Po(l - A) 



:i-a. 



*/(a+l) 



(12) 



a 3(a + l) 

and explicit expressions for p, p and V in terms of 0. Assuming a flat evolving universe 
described by the Friedmann equation H 2 = p (with H in units of H and p in units of 
Pent = 3Hq/8ttG), one obtains 

I . \y/l-A.(l 



3(a + l) 



In 



73(^+1) x 



where it is assumed that 
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v /l-A s (l-a^+ 1 )) + v ^X 
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One then readily finds the scalar field potential, 



In 



.i + v/t^X 



1/(0) = -Al +a p \ [cosh (3(a + 1)0/2)] ~ + [ cos h ( 3(a + 1)0/2)] ~ 



(13) 
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If one supposes that the energy density, p, may be decomposed into a mass varying 
CDM component, pm, and a dark energy component, p^, connected by the scalar field 
equations (jHl)-®, the equation of state (jBJ) is just assumed as an effective description of the 
cosmological background fluid of the universe. Since the CDM pressure, pm, is null, the 
dark energy component of pressure, p^, results in the GCG pressure, p = p^. Assuming 
that dark energy behaves like a cosmological constant, that is, its equation of state is given 
by p^{4>) = —P(j,{$), the dark energy density can be parameterized by a generic quintessence 
potential, p = U(<f>), since its kinetic component has to be null for a canonical formulation. 
It results in U(<f>) = —P(t>{4>) = P, where p is the GCG pressure given by Eq. (ITU1) . By 
substituting the result of Eq. ljTB"]) into the Eq. (fT0l) . and observing that H 2 = p, with p given 
by Eq.Q, it is possible to rewrite the GCG pressure, p, in terms of 0. It results in the 
following analytical expression for U(4>), 

'3(o + l) 



U{4>) = p^) = -p^) 



A„ cosh 



2a 
" 1+a 



(16) 



which is consistent with the result for V(4>) = (l/2)(p(0) — p{4>)) from Eq. ( 1X5]) . Since 
P<f>{<t>) + V^) — 0, the Eq. (jSJ) is thus reduced to 



dU(ct>) dp M dM(cb) 



d0 



dM d0 
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(17) 



and the problem is then reduced to finding a relation between the scalar potential U(<j>) and 
the variable mass M(4>). 

From the above equation, the effective potential governing the evolution of the scalar field 
is naturally decomposed into a sum of two terms, one arising from the original quintessence 
potential U(4>), and other from the dynamical mass M{4>). For appropriate choices of poten- 
tials and coupling functions satisfying Eq. (fl7|) . the competition between these terms leads 
to a minimum of the effective potential. In the adiabatic regime, the matter and the scalar 
field are tightly coupled together and evolve as one effective fluid. At our approach, once one 
assumes a A-like equation of state p^ = —p^, without any additional constraint on cosmon 
field equations, Eq. (1171) is naturally obtained. In the GCG cosmological scenario, the ef- 
fective fluid description is valid for the background cosmology and for linear perturbations. 
The equation of state of perturbations is the same as that of the background cosmology 
where all the effective results of the GCG paradigm are maintained. 

The Eq. ( fTBT) leads to p + p = pm + Pm which, in the CDM limit, gives 

p(a) + p(a) = M(a) n(a)(+p M EE 0). (18) 

Since the dependence of M on a is exclusively intermediated by <f){a), i. e. M(a) = M(4>(a)), 
from Eqs. (1101) and (1131) . after simple mathematical manipulations, one obtains 

2a 
1 + a 

(19) 

which is consistent with Eq. (fT7|) . From the above result, one can infers that the adequacy 
to the adiabatic regime is left to the mass varying mechanism which drives the cosmological 
evolution of the dark matter component. For mass varying CDM coupled with A-like dark 
energy, with p^ = —p^, the GCC leads to similar predictions for the equation of state, 
w = p/p, independently of the scale parameter a. The same is not true for hot dark matter 

(hdm) y. 

III. PERTURB ATIVE APPROACH FOR NEUTRINO MASSES 

Unlike photons and baryons, cosmological neutrinos have not been observed, so argu- 
ments about their contribution to the total energy density of the universe are necessarily 



M((j)) = M 



tanh (3(a + l)§) 
tanh (3(a + l)^) 
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theoretical. Otherwise, neutrinos coupled to dark energy can lead to a number of signifi- 
cant phenomenological consequences. The neutrino mass should be a dynamical quantity 
and neutrinos remain essentially massless until recent times. When their mass eventually 
increases close to its present value, their interaction with the background scalar field almost 

ceases 
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22]. The energy of the scalar field becomes the dominant contribution to the 



energy density of the Universe and cosmic acceleration ensues. 

Without loss of generality, such a behavior can be easily implemented to a degenerate 
fermion gas (DFG) of neutrinos. In the limit where T tends to in Eq. fl3]), the Fermi 
distribution f(q) becomes a step function that yields an elementary integral for the above 
equations, with the upper limit equal to the Fermi momentum here written as /3(a) = qp/a. 



It results in the equations for a DFG 



25], which can be useful for parameterizing the 



transition between ultra-relativistic (UR) and non-relativistic (NR) thermodynamic regimes. 
The equation of state can be expressed in terms of elementary functions of (3 = f3(a) and 
m = m(4>(a)), 



n(a) 



3tt 2 



Pm{a) = 



Praia) 



1 

8tt 2 

1 

87^ 



/3(2(3 2 + m 2 )^^ 2 + m? - arcsinh ((3/m) 



(20) 



)\/fi 2 + m 2 + arc sinh (/3/m) 



The relation given by Eq. (jl]) can be verified for the above definition. Besides the mass 
dependence, it is necessary to determine for which values of the scale factor the neutrino- 
scalar field coupling becomes important. The DFG approach take care of this. 

Let us then turn back to the explanation for the smallness of the neutrino masses. Ac- 
cording to the seesaw mechanism, the tiny masses, m, of the usual left-handed neutrinos are 
obtained through a very massive, M, sterile right-handed neutrino. The Lagrangian density 
that describes the simplest version of the seesaw mechanism via Yukawa coupling between 
a light scalar field and a single neutrino flavor is given by 



£ = m LR v L v R + M($)v R v R + h.c, 



(21) 



where it is shown that at scales well below the right-handed neutrino mass, one has the 
effective Lagrangian density [26] 



C 



m 



LR 



M{4>) 



VlV R + h.c. 



(22) 



Phenomenological consistency with the SM implies that logarithm corrections to the above 
terms are small, while it is well-know from the results of solar, atmospheric, reactor and ac- 
celerator neutrino oscillation experiments that neutrino masses given by m ~ m 2 LR (4>)/M((j)) 
lie in the sub-eV^ range. It is also clear that promoting the scalar field into a dynamical 
quantity leads to a mechanism in the context of which neutrino masses are time-dependent. 
Associating the scalar field to the dark e nerg y field allows linking NR neutrino energy den- 



sities to late cosmological times 



13, 



2l|. This scenario can be implemented through 



the perturbative approach via Eq. ((7j) [8]. It is evident that this approach is fairly general, 
as well as independent of the choice of the equation of state and of the dependence of the 
neutrino mass on the scalar field. 

The form of M(<j>) and of the equation of state can indeed lead to quite different scenarios. 
In particular, we consider three cases of neutrino mass dependence on 0. In case 01 we set 
m oc <f) 2 /M(4>), in case 02 we have set m oc 1/M(<£), and in case 03 we have set m oc (p~ 2 /M (</>). 
Energy densities, p, as function of the scale factor for the GCG fluid and for the decomposed 
components of the unified effective GCG background fluid, namely, the mass varying dark 
matter, p m , the cosmon-Zz&e dark energy, p^, and the perturbative DFG of neutrinos, p v , are 
obtained in the Fig. [TJ The neutrino densities are computed for neutrinos masses correlated 
to the dark matter mass, M{4>), in correspondence with the GCG scenario with A s = 3/4, 
that represents a realistic fraction of dark energy, and with a — 1, 1/2, 1/4 and 1/8. In 
the Fig. [2] we have taken into account the energy densities for photons (p 7 (a) = p 70 a" 4 ) 
and baryons (pb(o) = PBo a ~ 3 ) in complement to our effective GCG scenario. Then we have 
computed the relative modifications on the spectrum of the cosmic evolution of the dark 
sector components. 

Conveniently in agreement with the phenomenological predictions, we have set the present 
values of the energy densities as: = 0.95A 3 , p m0 = 0.95(1 — A s ), p u0 = 0.003, p 70 = 0.002, 

and pbo = 0.045. Assuming that neutrinos are NR at present, and observing the behavior 

i 

of Eqs. (12~TT) . we have set (3(a ) = (3(1) = pl . This condition allows one to establish the 
correspondence between the values of a for which the coupling transition takes place, i. e. 
the transition between relativistic and NR regimes, so that the neutrino masses achieve the 
present-day values. For our model such a transition is governed by the DFG behavior of 
Eq. (J2II). 

These scenarios are by no means the only possibilities. In particular, we choose them 
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as they correspond to the simplest feasible possibility for active neutrino mass generation, 
given that the scalar field has mass dimension one [26]. The essential information of the 
mass dependence on the scale factor is that, for case 03, the neutrino mass increases with 
decreasing 0, in opposition to what happens for cases 01 and 02. Through the stability 
analysis performed in the next section we shall elucidate some important consequences due 
to such a distinct behavior. 

In the Fig. [3] we verify how the energy density p and the corresponding equations of 
state w = p/p for the composed fluid deviate from the effective GCG scenario. For mass 
varying CDM coupled with dark energy with p^ = —p^, the effective GCC leads to similar 
predictions for w, independently of the scale factor a. The same is not true for HDM which, 
in the DFG approach, while weakly coupling with dark energy, leads to the same behavior 
of the GCG just at late times (a ~ 1). 



IV. ANALYSIS AND RESULTS 



The possibility of adiabatic instabilities in eosmo.og.ea. seenarios was previously pointed 
out [27] in a context of a mass varying neutrino model of dark energy. In opposition, in 
the usual treatment where dark matter are just coupled to dark energy, cosmic expansion 
together with the gravitational drag due to cold dark matter have a major impact on the 
stability of the cosmological background fluid. Usually, for a general fluid for which we know 
the equation of state, the dominant effect on the sound speed squared c 2 s arises from the 
dark sector component and not from the neutrino component. 

For the models where the adiabatic regime (cf. Eq. ( !T7|) ) implies an equation of state 
reproducing cosmological constant, A, effects, p^ = — p^, one obtains c 2 s = —1 from the 
very start of the analysis. The effective GCG is free from this inconsistency. A detailed 
quantitative analysis of the stability conditions for the GCG in terms of the squared speed 
of sound is discussed in Ref . [l3| . Positive c 2 implies that < a < 1 in the GCG equation of 
state. The coupling of the dark energy component with dynamical dark matter is responsible 
for removing such inconsistency by setting c 2 s ~ ^ > 0. 

In the Fig. H] we show the results for c 2 for a cosmological background scenario corre- 
sponding to the sum of the energy components of an effective GCG fluid: mass varying dark 
matter (DM), cosmon-Zz/ce dark energy (DE) and perturbative DFG of neutrinos (u). We 
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have considered the GCG phenomenological given by A a = 3/4 and a = 1 (solid line), 1/2 
(dot line), 1/4 (dash dot line), and 1/8 (dash line). The perturbative influence of neutrinos 
on the positiveness of c 2 s has its magnitude measured by comparing the results with those 
for the GCG scenario with and and without neutrino perturbation for each one of the three 
cases in correspondence with Fig. [TJ In this scenario, the most important result illustrated 
by Fig. H] is that growing neutrino masses parameterized by <p 2 /M(4>) (case 03) results in 
instabilities at present times. 

This instability is characterized by a negative squared speed of sound for an effective 
coupled neutrino/dark energy fluid, and results in the exponential growth of small scale 
modes 6|. A natural interpretation for this is that the Universe becomes inhomogeneous 
with the neutrinos forming denser structures or lumps jl^l H . 29]. Effectively, the scalar 
field could mediate an attractive force between neutrinos leading to the formation of neutrino 
nuggets. Anyway, these results are consistent with the accelerated expansion of the universe 
ruled by the dynamical masses of CDM and neutrinos since positive values for (1 + 3(p^ + 
Vm)/{.P<j> + Pm)) are observed, as we can notice in the Fig. El 



V. CONCLUSIONS 



Growing neutrinos coupled to mass varying dark matter plus cosmon- Me dark energy 
were studied assuming that the cosmological background unified fluid presents an effective 
behavior similar to that of the GCG. The essential ingredient of this class of models is a 
neutrino mass that, through the seesaw mechanism, depends on the cosmon field and grows 
in the course of the cosmological evolution. Our setup is an effective GCG decomposed into 
two interacting components |8j. The first one behaves like mass varying CDM since it is 
pressure- less. The second one corresponds to a cosmon- like dark energy component <fi with 
equation of state given by p^{4>) = — p<f>{4>)- Apparently the model does not look different 
from the interacting quintessence models where one has two different interacting fluids. The 
present fraction of dark energy is set by a dynamical mechanism. As soon as dark matter 
become NR, their coupling to the cosmon triggers an effective stop (or substantial slowing) 
of the evolution of the cosmon. Before this event, the quintessence field follows the tracking 
of the minimum of an effective potential, in a kind of stationary condition which drives the 
adiabatic regime, and for which the coupling strength is strong compared to gravitational 
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strength. 

In the scope of finding a natural explanation for the cosmic acceleration and the cor- 
responding adequation to stability conditions, we have proposed a systematic procedure 
to treat variations of fundamental parameters of neutrino cosmology (0, m(4>) and M(<p)) 
independent of any particular theoretical model which enforces restrictive relations among 
these parameters. We have quantified the perturbative influence of such variations on the 
calculation of the squared speed of sound and on the cosmic acceleration of the universe, 
noticing that the former quantity is much more sensible to variations with respect to the 
exact GCG background cosmology. 

An interesting feature is however that the GCG has its cosmological evolution reproduced 
by a CDM mass varying mechanism consistent with the mass generation mechanism for 
neutrinos. Most of dark sector models predict an entangled mixture of interacting dark 
matter and dark energy. We have consistently introduced the possibility of such a mixed 
coupling with active and sterile neutrinos, for which we have tested three hypothesis for 
seesaw masses. 

At our approach, the mass scale M and scalar field are tightly coupled together and 
evolve as a unique effective fluid. The effective potential driving the evolution of the scalar 
field is decomposed into a sum of two terms, one arising from the original quintessence po- 
tential U(4>), and the other from the dynamical mass of the dark matter. The sterile neutrino 
mass produces the mass varying dark matter effects which indirectly result in MaVaN's, and 
the active neutrino has its mass perturbatively computed via seesaw relations. Since the 
results for the cosmic acceleration followed by cosmological stability for GCG scenarios are 
well-defined, one can assert that the increase of neutrino mass acts as a cosmological clock 
or trigger for the crossover to the effective scenario here studied. In fact, a scalar field asso- 



ciated to dark energy in connection with the SM neutrinos and t 
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re electroweak interactions 



321 ]. The case of cosmolog- 



may bring important insights on the physics beyond the SM 
ical neutrinos, in particular, is a fascinating example where salient questions concerning SM 
particle phenomenology can be addressed and hopefully better understood. 
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FIG. 1: Energy densities as function of the scale factor for the GCG fluid (GCG) and for the 
decomposed components of the unified background fluid (effective GCG): mass varying dark matter 
(DM), pm, cosmon- like dark energy (DE), p^^ and perturbative DFG of neutrinos, Pm- 

For the 03 

cases of MaVaN's, the dynamical mass dependencies on <fi are: 4> 2 /M ((/>), 1/M((f>) and <p~ 2 /M ((/>). 
Neutrino densities are computed for a mass varying behavior of M{4>) in correspondence to a GCG 
scenario with A s = 3/4 and a = 1 (solid line), 1/2 (dot line), 1/4 (dash dot line), and 1/8 (dash 
line). 
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FIG. 2: Energy densities (p/p Clit ) as function of the scale factor for all energy components of the 
universe. The results are compared with those of an exact GCG scenario with A a = 3/4 and 
parameters a = 1 (solid line), 1/2 (dot line), 1/4 (dash dot line), and 1/8 (dash line). We have 
considered all the 03 cases for MaVaN's in correspondence with Fig. [U 
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FIG. 3: Equation of state w = p/p for 03 cosmological background scenarios in dependence on a. 
Notice the separate behavior of the neutrino fluid (magenta) and its influence on the equation of 
state for the GCG scenario Equations of state for decoupled GCG (black) and for GCG coupled 
to neutrinos (blue) are compared for each one of the 03 cases in correspondence with Fig. [TJ The 
perturbative influence of neutrinos on the GCG background scenario is negligible. The results are 
for a GCG scenario with A s = 3/4 and parameters a = 1 (solid line), 1/2 (dot line), 1/4 (dash dot 
line), and 1/8 (dash line). 
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FIG. 4: Squared speed of sound c 2 s = dp/dp for 03 cosmological background scenarios in dependence 
on a. Notice the results for the decoupled neutrino fluid (magenta). The perturbative influence of 
neutrinos on the positiveness of c 2 s has its magnitude measured by comparing the curves for the 
GCG scenario with (blue) and and without (black) neutrino couplings for each one of the 03 cases 
in correspondence with Fig. [TJ Even for unstable neutrino scenarios (third plot), the perturbative 
influence on the GCG background scenario is negligible. The results are for a GCG scenario with 
A s = 3/4 and parameters a = 1 (solid line), 1/2 (dot line), 1/4 (dash dot line), and 1/8 (dash 
line) . 
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FIG. 5: Cosmic acceleration rate a/{ap) for 03 cosmological background scenarios in dependence 
on a. We have considered the GCG (black), the GCG plus neutrinos (blue), and the GCG plus 
neutrinos, baryons and radiation (red). Notice that the cosmic acceleration is essentially driven by 
the dark sector corresponding to the GCG. The perturbative influence of neutrinos on the cosmic 
acceleration is not relevant. We use the same GCG parameters from Fig. [U 
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